Intercalant and intermolecular phonon assisted superconductivity in K-doped picene 
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K3-picene is a superconducting molecular crystal with critical temperature T c = 7K or 18K, 
depending on preparation conditions. Using density functional theory we show that electron-phonon 
interaction accounts for T c 3 — 8 K. The average electron-phonon coupling, calculated by including 



the phonon energy scale in the electron-phonon scattering, is A = 0.73 and 0J\ o 



18.0 meV. 



Intercalant and intermolecular phonon-modes contribute substantially (40%) to A as also shown by 
the isotope exponents of potassium (0.19) and carbon (0.31). The relevance of these modes makes 
superconductivity in K-doped picene peculiar and different from that of fullerenes. 

PACS numbers: 31.15.A-, 74.70.Kn, 63.20.kd 
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Most materials made of organic molecules are not 
metals as the valence bands are completely filled and 
separated by the conduction band. However, the dis- 
covery in 1980 of superconducting tetra-methyl-tetra- 
selenium-fulvalene[l[, demonstrated the possibility to in- 
duce superconductivity by intercalation of large organic 
molecules with inorganic anions, serving as acceptors or 
donors. This general rule was confirmed by the discov- 
ery of superconductivity in alkali-doped fullerenes 0] with 
transition temperatures (T c ) up to 40 K in CS3C60 0- 

The field was renewed recently with the discovery of su- 
perconductivity in potassium (K)-intercalated molecular- 
solids based on polycyclic aromatic hydrocarbons, such 
as picene 0] (T c = 7K or T c = 18K Q), phenantreneQ 
(T c = 5K) and coronene (T c = 15K). These systems 
are very appealing as they exist in a large variety, and it 
is possible to tune their properties in many ways. For in- 
stance, the chemistry and functionalization of polycyclic 
hydrocarbons is well known and relevant for environmen- 
tal and medical issues Q, the intercalation with alkali 
atoms is often possible, and the electrochemical doping 
in picene field-effect transistors has been achieved [ij]. 

In addition, organic molecular crystals are particu- 
larly appealing for technological applications such as or- 
ganic light emitting diodes (OLED), field effect transis- 
tors (OFET) or solar cells. Consequently, the study of 
charge injection and carrier transport in these materi- 
als are central to the understanding and improvement 
of such devices. In particular, the study of the electron- 
phonon coupling interaction is crucial to describe the dif- 
ferent regimes of charge transport [loj . 

In molecular crystals, the electrons can couple to three 
kinds of vibrations; intramolecular, intermolecular (rigid 
translations and rotations of molecules) and intercalant 
jjhonons. In alkali-doped Cgo, after a long debate, |3.HlT- 
14| . it has been established that superconductivity is 
driven mainly by intramolecular modes. This conclusion 
is also supported by the negligible Rb isotope effect mea- 
sured in RbaCeo HI However the situation could be 



substantially different in K-doped Picene where a more 
prominent role of intercalant modes in superconductivity 
is suggested by the significant hybridization of the low- 
energy bands with K atoms observed in density func- 
tional theory (DFT)[l3|. Indeed, the electronic struc- 
ture of K-doped picene is different from that obtained 
by a rigid-doping of the picene crystal without explicit K 



ionslT7j. Thus, at the moment, it is still unclear (i) if su- 
perconductivity in K-doped picene is phonon-mediated, 
(ii) if the intramolecular phonons alone provide enough 
coupling or if intercalant and intermolecular vibrations 
should be also considered. 

These questions were first addressed in an early DFT 
calculation where it was found a sizable electron- 
phonon coupling, dominated by high-energy intramolec- 
ular in-plane carbon (C) vibrations at 1000 cm -1 . How- 
ever, this work uses a rigid-band doping approximation 
of the pristine picene crystal structure and thus neglects 
the molecular relaxation after intercalation, the K contri- 
bution to the electron-phonon coupling, and the screen- 
ing of the electron-phonon deformation-potential by the 
metallic bands. Furthermore only modes above 100 cm -1 
were considered. 

In this letter, we perform first-principles DFT calcu- 
lations on the superconducting properties of K3-picene. 
Removing all possible a-priori approximations, and in 
contrast to [l8j, we demonstrate the critical role played 
by the intercalant and intermolecular modes in the 
electron-phonon coupling mechanism to account for the 
experimental T c . We give a numerical prediction of the 
isotope exponents in the molecular crystal and relate 
them to the phonon-mediated superconductivity. 

We compute the electronic structure of K3-picene 
within DFT in the local density approximation 
(LDA)fl9|. We relax the molecular geometr y in the P2i 
symmetry in the experimental unit cell 0, 21 1, obtain- 
ing a structure similar to that of [171 ] with the molecules 
arranged in a herringbone configuration and the inter- 
calants aligned parallel to the molecular planes [22j|. 
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FIG. 1: DFT-LDA band structure of K3-picene near the 
Fermi level. The DOS per spin per cell (comprising 2 picene 
molecules) is shown in the right panel. Note the strong DOS 
variation in a range of ±hoji og from the Fermi energy, with 
the characteristic phonon frequency uji og = 18 meV 



The band structure of the K3 picene is reported in 
Fig. [TJ The molecular nature of the compound gives rise 
to very narrow bands ek n , with bandwidth Ae of few 
tenths of eV. The Fermi level e F in the K-doped picene 
crosses the bands originating from the LUMO+1 (sec- 
ond Lowest Unoccupied Molecular Orbital) states of the 
undoped molecule, which are entangled to those com- 
ing from the LUMO+2 states, and very close to the 
ones belonging to the LUMO bands. The bandwidth 
Ae has the same order of magnitude as the characteristic 
phonon frequency uj p h [the logarithmic phonon average 
Wi og =18 meV, calculated in this work (see below)]. The 
density of states (DOS) varies substantially in the range 
[e F - ioiog, e F + wiog], as shown in Fig. [TJ 

We carry out phonon calculations in the density func- 
tional perturbation theory framework (2^ . We consider 
phonon momenta on a JV g = 2x2x2 grid. We get 
dynamically stable phonons, signaling a well converged 
phonon dispersion w q „ and no structural instabilities in 
the chosen geometry. In FiglJJc) we report the projected 
phonon DOS, p s (e) = J] qj/ e qj/ • Vse^luj^ - e), where 
the phonon eigenstates e q „ are 3N-dimensional vectors, 
with N the number of atoms, and the projector Vs is a 
3Nx3N tensor. We split the full phonon space into potas- 
sium (K), hydrogen (H), in-plane (C = ) and out-of-plane 
Cx carbon S'-subspaces (the plane is the one contain- 
ing each picene molecule). We also consider partition- 
ing into K, intermolecular and intramolecular vibrations. 
The projections reveal the presence of the K and inter- 
molecular modes at low frequencies, while only a very 
small fraction of intramolecular modes is present below 
200 cm -1 . On the other hand, the in-plane C and H 
modes occupy the high-energy part of the DOS. We com- 
pute the phonon spectrum also for the neutral isolated 
molecule, Figj2jd) . The low-energy spectrum of molecule 
and crystal differ considerably. In the solid the rigid- 
body rototranslations (intermolecular modes) are spread 
up to 300 cm^ 1 , while in the molecule they are pinned 



at zero frequency. In addition we observe in the crystal 
a significant spectral weight transfer of the low-energy 
out-of-plane C oscillations to higher frequencies. 

0.80 F ' ' ' ' ' ' 1 <-^\ 




200 400 600 800 1000 1200 1400 1600 

co (cm -1 ) 



FIG. 2: Panels (c) and (d): phonon DOS of the Fv3-picene 
crystal and picene molecule, respectively. The DOS is pro- 
jected on phonon subspaces S which correspond to hydrogen 
(H), potassium (K), out-of-plane and in-plane carbon (Cx 
and C=) modes. The black solid curve is the sum of K and 
intermolecular modes. In the (d) panel the black line is absent 
as there is no K and all phonons are intramolecular, except 
for the six rototranslational modes at zero frequency. We do 
not show the C-H stretching modes, which do not couple to 
electrons, and are located around 3050 cm -1 , well above the 
frequency range plotted. Panels (a) and (b): X(u) anda 2 F(u>) 
resolved by projections on H, K, Cx, and C= phonons. A(cj) 
ramps up by 90% in the first 300 cm -1 . The black solid line 
is the sum of K and intermolecular modes, which contribute 
to 40% of the total A. 

For each phonon mode v with momentum q we com- 
pute the electron-phonon interaction: 

= W 2 JV(0)JV fc ^ ISkn,k+ q J 2 
4 7 k,n,m 

X (jfkn — ik+q,m) <5(ek+q,m — 6kn — Wqi/), (1) 
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that couples the occupied state |k, n) of momentum k and 
band n with the empty state |k + q, m) separated by the 
phonon energy w q „. The k-summation in Eq. Q]is over 
the Brillouin zone on a N k electron-momentum mesh, 
and N(0) (= 6.2 eV _1 ) is the electron DOS per spin per 
cell at the Fermi level. The electron-phonon matrix ele- 
ments are 5k„,k+ q ™ = E s (k,k + q)/ A /2M s cj qi/ , 
where e q!/ is the 3-dimensional c\y eigenphonon compo- 
nent on the s-th atom with mass M s , and d* m (k, k+q) = 
(k + q,m\SV/Suq S \k,n), with 5V/Suq S the deformation 
potential for a given phonon displacement qs of the s- 
th atom. In Eq. [T] /kn are Fermi functions depending 
on the temperature T, and the expression for A qj/ has 
been evaluated by a T — »■ extrapolation. In the "adi- 
abatic" limit, namely for <C Ae, where Ae is the 
bandwith, the expression for A q „ in Eq. [1] reduces to the 
one proposed by Allen [5^], and generally used in previous 
electron-phonon estimates: 



\AD _ f 

uj^N(0)N k 



E 

k.n.m 



^kn,k+qm| 2 ^( e k,n)5(ek+ qi m)- 



(2) 

We are going to dub X AD in the Equation above as "adi- 
abatic" , while A in Eq. [T] as "non-adiabatic" . 

The electron-phonon matrix elements g£ n k+qm are 
computed on a highly dense electron-momentum grid 
by means of the recently developed Wannier interpola- 
tion formalism [26M30j ] . This allows one to perform con- 
verged calculations with a delta function Gaussian smear- 
ing a much smaller than the bandwidth Ae fsilj. Sur- 
prisingly, the exact zero-smearing extrapolation leads to 
a strong electron and phonon-momentum dependence 
in g, and so in A and X AD . This is in contrast 
with the conventional molecular crystal picture, where 
g and the electron-phonon coupling are independent of 
the electron and phonon momenta. To quantify this 
dependence, we evaluated the coupling strength aver- 
age per molecule, given by 2A q £) A r (0)/J q , where J q = 
E k , n ,m S(e K n)S{ew+q, m )/N k is the nesting factor, \ AD is 
summed over the phonon modes, and the factor of 2 ac- 
counts for the number of molecules in the unit cell. One 
obtains the phonon-momentum integrated average- value 
of 172 meV, with large variations of ±40% as a function 
of the phonon-momentum [23|]. A detailed analysis of 
the dependence of A on q and er is presented in [23j. 

In Tab.U we report our results for the electron-phonon 
coupling. There is a significant difference between the 
"adiabatic" and "non-adiabatic" values, as uj p h ~ Ae. 
Therefore, one must keep the full expression in Eq. [T]for 
an accurate estimate of A. In Tab. [H we also report the 
phonon frequency logarithmic average (oJiog)- 

An estimate of T c is given by the McMillan for- 
mula [32], which requires A, un og , and the screened 
Coulomb pseudopotential fi* . With /i* G [0.1,0.2], we 
obtain T c in the range 3 — 8 K, which includes the value 
of 7 K reported in [J] for one of the two phases of K3- 





A 


wiog (meV) 


T C (K) 
fi* 6(0.1,0.2] 


«(C) 


a(K) 


a(H) 


adiabatic 


0.88 


25 


8-16 


0.32 


0.13 


0.05 


non-adiabatic 


0.73 


18 


3-8 


0.31 


0.19 


0.00 



TABLE I: A, phonon frequency logarithmic average uii og , 
McMillan critical temperature T c (obtained with in the 
range of 0.1-0.2), and isotope exponents (for fi* — 0.14), com- 
puted based on the electron-phonon coupling expression in 
Eq. [TJ (non-adiabatic), which explicitly includes the phonon 
energy ui in the electron-phonon scattering process, and the 
one in Eq. [2] (adiabatic), which is a Wpd/Ae — > approxi- 
mation of Eq. [1] The "non-adiabatic" values are our best 
estimates of the various quantities. 



picene [5j . Caution must be taken in using the McMillan 
formula, as the on-site molecular correlation can be im- 
portant with respect to the bandwidth[18, 33], and also 
"non-adiabatic" effects in the vertex corrections can in- 
validate the Migdal's theorem[34|. However, our estimate 
of T c is an indication that the electron-phonon coupling is 
strong enough to explain the mechanism behind the man- 
ifestation of the superconducting phase. We computed 
also the isotope exponents a(X) = —d\ogT c /d\ogM x 
(see Tab. Q} for all constituents of the molecular crys- 
tal. It turns out that the role of H is negligible, while 
the large exponent of potassium (K) points toward the 
important role played by the intercalant in the electron- 
phonon coupling to set the value of T c . 

For a given value of /x* our predicted T c are more than 
8 times smaller than that predicted with A and ui og found 
in [lij]. This comes mainly from the value of 0Ji og (1021 
cm -1 ) in fisj j . very different from our best estimate of 18 
meV (145 cm -1 ). This large discrepancy may come from 
the intercalation driven structural change not considered 
in that work. Indeed, the electronic DOS of Fv3-picene, 
computed in the rigid-doping approximation with the un- 
doped geometry, is 14.12 eV _1 per spin per cell frR flij ]. 
a value 2.3 larger than that found here and in [17j for 
the relaxed geometry with K. In addition, [l8[ neglects 
the screening of the electron-phonon interaction by the 
metallic bands, the contribution of K phonons and vibra- 
tions of energy lower than 100 cm . 

To understand our findings, we compute the Eliash- 
berg function a 2 F(uj) = Y,q V A q „w q „ 8(co - uj clv )/(2N q ), 
plotted in Fig. [2][b), and the integral X(uj) = 
2 J Q duj' a 2 F '(a/) / uj' , shown in Fig. [2ja), which gives the 
total electron phonon coupling A = ^q^/^q m the 
lu — > 00 limit. We note that A(w) converges very rapidly 
to a large fraction of the total value in the first 300 cm -1 . 

We then decompose A into a set of phonon subspaces 
S, as we did for the phonon DOS, as in [35j. The evalua- 
tion of the projected electron-phonon matrix elements 

9s = E s (PsV) s ' <m( k > k + q)/^ 2M ^q- follow s 
straightforwardly, and the corresponding projected Azjf 
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reads: 



*qz; 



w\ v N(0)N k 93 9s ' 

H k,n,m 
x (/kn — /k+q,m) <S(ek+q,m — £kn — ^qv)) (3) 



where for the sake of readability we dropped out the mo- 
mentum and orbital indexes from g$. A = ff S ~ 
J2s S' Sq v > an£ i the contribution of each subspace 
S is computed as ^2 S , X S ' S , where we add both the di- 
agonal term and the usually very small off-diagonal con- 
tributions. The sum over K and intermolecular modes 
is reported in Fig. [2] as a black solid line. It gives more 
than 40% of the total A, with half of it due to the alkali 
intercalant, and the corresponding a 2 F{u>) localized at 
small frequencies (up to 300 cm -1 ). The intramolecu- 
lar phonons contribute to the remaining 60%, with the 
main role played by the out-of-plane C modes. This 
is a remarkable result, and shows that K-doped picene 
is far from being a prototype molecular crystal. The 
electron-phonon coupling is only partially supported by 
intramolecular phonons, and the q-dependence of the 
electron-phonon coupling is strong. 

To summarize, we found that in the K3-picene the 
electron-phonon coupling is large enough to explain the 
experimentally measured T c of 7 K. Despite the molec- 
ular nature of the crystal, a significant contribution to 
the coupling is given by the dopant and intermolecu- 
lar phonon modes, which account for 40% of the total 
A. This has a strong impact on the isotope exponent 
of potassium, whose value turns out to be large (close 
to 0.20), and represents an experimentally-accessible sig- 
nature of the importance of non-intramolecular modes. 
Our work shows that there are fundamental differences 
between the families of picene and fullerene. The demon- 
stration of the importance of coupling between the elec- 
trons at the Fermi level with intermolecular and inter- 
calant low-energy phonon modes will have also a rel- 
evant contribution to the developing of transport the- 
ories, usually accounting for the local electron-phonon 
coupling (intramolecular coupling) and considered in a 
small-polaron theory of narrow bands. Extensions to in- 
clude non-local couplings (intermolecular coupling) are 
required [1^]. This opens the way to a more compre- 
hensive view on superconductivity, transport, and other 
phenomena in the larger and larger group of recently dis- 
covered materials based on metal-intercalated aromatic 
hydrocarbons. 

We acknowledge Y. Kubozono for useful discussions. 
Computer allocation has been supported by CINECA- 
HPC ISCRA and EU DEIS A-SUPERM AG . 
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GEOMETRY 

We report the coordinates of the 78 atoms in the monoclinic unit cell with a,b,c = 8.707, 5.912, 12.97A, (3 = 92.77°. 
In the Table, they are split into two groups, one for each K 3 -picene molecule, related each other by the P2i symmetry 
operations. The atomic symbols are in the first column. The coordinates are in crystal units. 





molecule 1 


molecule 2 


c 


0.4556221 





6985366 





5225742 


-0 


4556221 


0.1985366 


-0 


5225742 


c 


0.3666947 





7916786 





5978784 


-0 


3666947 


0.2916786 


-0 


5978784 


c 


0.4078295 





6965432 





4186143 


-0 


4078295 


0.1965432 


-0 


4186143 


H 


0.5603062 





6067382 





5464734 


-0 


5603062 


0.1067382 


-0 


5464734 


c 


0.2680166 





7885757 





3834603 


-0 


2680166 


0.2885757 


-0 


3834603 


H 


0.4802802 





6105565 





3645861 


-0 


4802802 


0.1105565 


-0 


3645861 


C 


0.1758863 





9015619 





4566338 


-0 


1758863 


0.4015619 


-0 


4566338 


c 


0.2192085 





7895931 





2754507 


-0 


2192085 


0.2895931 


-0 


2754507 


c 


0.0759228 





8943093 





2435217 


-0 


0759228 


0.3943093 


-0 


2435217 


c 


0.3100630 





7014734 





1996684 


-0 


3100630 


0.2014734 


-0 


1996684 


c 


0.2617866 





6976639 





0950588 


-0 


2617866 


0.1976639 


-0 


0950588 


H 


0.4193296 





6208887 





2225564 


-0 


4193296 


0.1208887 


-0 


2225564 


c 


0.1210633 





7886946 





0646123 


-0 


1210633 


0.2886946 


-0 


0646123 


H 


0.3366770 





6247067 





0390565 


-0 


3366770 


0.1247067 


-0 


0390565 


C 


0.0298681 





8877619 





1371182 


-0 


0298681 


0.3877619 


-0 


1371182 


H 


0.0819831 





7857442 


-0.0168198 


-0 


0819831 


0.2857442 


0.0168198 


H 


-0.0771974 





9716713 





1101440 


0.0771974 


0.4716713 


-0 


1101440 


C 


-0.0065371 


1 


0089684 





3165922 


0.0065371 


0.5089684 


-0 


3165922 


C 


0.0469952 


1 


0175338 





4200294 


-0.0469952 


0.5175338 


-0 


4200294 


H 


-0.1098166 


1 


1064455 





2929669 


0.1098166 


0.6064455 


-0 


2929669 


H 


-0.0159085 


1 


1243530 





4725799 


0.0159085 


0.6243530 


-0 


4725799 


C 


0.2237032 





8990497 





5655413 


-0 


2237032 


0.3990497 


-0 


5655413 


C 


0.1396179 


1 


0059855 





6385974 


-0 


1396179 


0.5059855 


-0 


6385974 


c 


0.1916193 


1 


0169922 





7421052 


-0 


1916193 


0.5169922 


-0 


7421052 


H 


0.0315715 


1 


0922925 
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BROADENING AND q-DEPENDENCE OF A 

In the Figure, we present the phonon momentum q-resolved A's, evaluated as in Eq. 2 in the paper, as a function 
of the Gaussian smearing a of the delta functions. The plotted results are obtained by using a 120 3 electron k-point 
mesh. The calculations performed with 60 3 k-point mesh give identical results down to T =150K and a =4.3 meV, 
with T — 3c The q vectors are reported in crystal units. The exact zero-smearing extrapolation leads to a strong 
q-dependence in A, which can be seen only for a much smaller than the bandwidth Ae (<C 0.1 eV). In the inset, we 
show the q-integrated "adiabatic" and "non-adiabatic" A. The "non-adiabatic" A is extrapolated with respect to both 
a and T. 

Note the divergence of at the zone center. This occurs in general in presence of optical modes, because of the 
divergence of the nesting factor J q at q=0. For this reason, in the "adiabatic" q-averaged X AD the T point has been 
excluded. In the "non-adiabatic" case, A q is well behaved also at T. 
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COUPLING STRENGTH AVERAGE 

From Aq 15 , we can quantify the q-dependent coupling strength average per molecule, by the formula 2A q ' D -/V(0)/J q , 
where J q = ^ k n m 5(ek,n)^(ek+q.m)/-A r fc is the nesting factor, A q D is summed over the phonon modes, and the factor 
of 2 accounts for the number of molecules in the unit cell. One obtains the values (in meV) { 183, 193, 131, 133, 239, 
146, 191, 161 } for the q-vectors (in crystal units) { (0,0,0), (0,0, i), (0,|,0), (0,±,±), (5, 0,0), (± 0, ±), (3,5,0), 
i> 5) I' respectively. Note that even at V the ratio is well defined, as the divergences of both Ap and Jr cancel 
out. 



